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Abstract 

Using the information theory approach, in both its extensive and nonextensive versions, we es- 
timate the inelasticity parameter K of hadronic reactions together with its distribution and energy 
dependence from pp and pp data. We find that the inelasticity remains essentially constant in energy 
except for a variation around K ~ 0.5, as was originally expected. 

PACS numbers: 96.40 De, 13.85.Tp, 96.40.Pg, 13.85.Ni 

1 Introduction 

The inelasticity K (0 < K < 1) of a reaction has well established importance in working with data from 
cosmic ray cascades pQ (cf. also El E] and references therein) . It tells us what fraction of the energy of 
a projectile is used for production of secondaries and what fraction flows further along the cascade chain. 
In cosmic ray observables, K in fact appears in some combination which also contains the mean free path 
for the particle propagation in the atmosphere, or equivalently, the total inelastic cross section Ui n . This 
makes K difficult to estimate because of the freedom available to attribute the observed effects either 
to K(s) or to cr(s). It was therefore proposed in J that in order to extract K(s) unambiguously from 
cosmic ray experiments one should analyse simultaneously the data from (at least) two different types of 
experiments for which combinations of K and cr;„ arc different 5 . 

Nowadays there is a strong tendency to replace inelasticity and simple energy flow models by more 
refined and complicated models of multiparticle production (see [3111101 and references therein). In their 
present formulations such models differ substantially among themselves, concerning both their physical 
basis and the (usually very large) number of parameters used, and lead to quite different, sometimes even 
contradictory, predictions |S] . While developing models is necessary for the global understanding of 
cosmic ray physics, for the purpose of studying energy flow it may be desirable to have a more economical 
description of high energy collisions, involving only a small number of parameters. This is one of the 
advantages of working with the concept of inelasticity ^ ^] . 

The inelasticity is also a very important quantity in phenomenological descriptions of hadronic and 
nuclear collisions in terms of statistical models of multiparticle production processes |11II12| . In this case 
it enters either explicitly, as a single parameter K defining the initial energy M — K-^fs to be hadronized, 
or implicitly, when one finds that out of numerous parameters of the model only the combination leading 
to the given fraction of available energy to be converted into produced secondaries, K ■ ^fs, is important, 
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see Attempts to estimate it are thus fully justifiable. In ^2 113 the calculations of the inelasticity 
distribution, x(K) (and its energy dependence), based on the assumed dominance of high energy multi- 
particle processes by gluonic interactions, was presented. In the other calculations the mean inelasticity 
K and its possible energy dependence was simply estimated either by using thermal-like model formulas 
applied to collider pp data |14l 1151 ITS1 1171 118) (like, for example, in J^j) or by some other means |2(J) . 

In this paper we address this problem again, this time by means of the information theory approach 
both in its extensive [2U and nonextensive [221 versions. The idea is to describe the available data 
by using only a truly minimal amount of information avoiding therefore any unfounded and unnecessary 
assumptions. This is done by attributing to the measured distributions (written in terms of the suitable 
probability distributions) an information entropy S and maximizing it subject to constraints which ac- 
count for our a priori knowledge of the process under consideration. As a result one gets the most probable 
and least biased distribution describing these data, which is not influenced by anything else besides the 
available information. In such approach the inelasticity K emerges as the only real parameter; all other 
quantities being well defined functions of it. We attempt to clarify here the role of the inelasticity by 
using both the extensive and nonextensive versions of information theory. All necessary background on 
the information theory approach needed in the present context is given in the next Section. Section 3 
contains our results for the pp and pp collisions. Our conclusions and summary are presented in the last 
Section. 



2 General ideas of information theory approach 

As presented at length in |21j (where further details and references can be found) the information theory 
approach provides us, by definition, with the most probable, least biased estimation of a probability 
distribution {pi, i = 1, . . . , n} using only knowledge of a finite number of observables {Fk, k = 1, . . . , n} 
of some physical quantities obtained by means of pi and defined as: 

n 

F k = (F k )=J2p l Ft t) - (1) 

i=l 

One is looking for such {p{\ which contain only information provided by {F^} and nothing more, i.e., 
which contain minimal information. The information connected with {pi} is quantified by the Shannon 
information entropy defined as: 

S = -Ej^lnpi. (2) 

Minimum information corresponds to maximum entropy S, therefore the {pi} we are looking for are 
obtained by maximizing the information entropy S under conditions imposed by the measured observables 
{Fk} as given by eq.Q). They result in a set of Lagrange multipliers {A*,, k — 1, . . . , r} and the generic 
form of {pi} we are looking for is |21| : 

(3) 

where Z is obtained from the normalization condition 53i=iP« = 1- 

Such approach was applied long time ago to experimental data on multiparticle production with the 
aim at establishing the minimum amount of information needed to describe them |13| . The rationale 
was to understand what makes all the apparently disparate (if not outright contradictory) models of that 
period fit (equally well) the data. The result was striking and very instructive |13j : the data considered 
(multiplicity and momentum distributions) contained only very limited amount of information, which 
could be expressed in the form of the following two observations: (i) the available phase space in which 
particles are produced is limited (i.e., there is some sort of pr cut-off) and (ii) only a part K £ (0, 1) of 
the available energy y/s is used to produce the observed secondaries, the rest being taken away by the so 
called leading particles (i.e., inelasticity emerges as one of the cornerstone characteristics of reaction). All 
other assumptions, different for different models (based on different, sometimes even contradictory, phys- 
ical pictures of the collision process) were therefore spurious and as such they could be safely dropped out 
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without spoiling the agreement with experimental data. In fact, closer scrutiny of these models showed 
that they all contained, explicitly or implicitly, precisely those two assumptions mentioned above and 
that was the true reason of their agreement with data. 

In this paper we are therefore following the same line of approach with the aim at deducing from 
the available data for the inelasticity parameter K. Notice that the formula resembles the statistical 
model formulas based on the Boltzmann-Gibbs statistics as used in However, in J3J) no thermal equi- 
librium is assumed and all Afc (i.e., among others also the "partition temperature" T in are given by 
the corresponding constraint equation JQ) whereas normalization fixes Z, which is a free parameter in |19|. 

This approach can be generalized to systems which cannot be described by Boltzmann-Gibbs (BG) 
statistics because of either the existence of some sort of long-range correlations (or memory) effects or 
fractal structure of their phase space or because of the existence of some intrinsic fluctuations in the 
system under consideration. It turns out that such systems are nonextensive and therefore must be 
described by a nonextensive generalization of BG statistics, for example by the so called Tsallis statistics 
|2"2") defined by the following form of the entropy: 

S q = - j-^SiCl-p?). (4) 
It is characterized by the nonextensivity parameter q such that for two independent systems A and B 

S q (A+B) — S q A + S q B + (1 — q)S q AS q B- (5) 

Notice that in the limit q — > 1 one recovers the previous form of Boltzmann-Gibbs-Shannon entropy J5J. 
Maximazing S q under constraints, which are now given in the form |24j : 

n 

Fk 9) = {F k ) q = J2\Pi] q Fk ,i \ (6) 
i=i 

results in the following power-like form of the (most probable, least biased) probability distribution: 



(<?) 1 
Pi = Pi = -77- exp 

Z 9 



where Z q is obtained from the normalization condition y^,-_ 1 p\ — !• and where 
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(7) 
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Of special interest to us here will be the fact that intrinsic fluctuations in the system, represented by 
fluctuations in the parameter 1/A in the exponential distribution of the form ~ exp(— x/A) result in its 
nonextensivity with parameter q given by a normalized variation of fluctuation of the parameter 1 /A 23 : 



q= 1± X>A; / n — ■ (9) 

<i> 2 

So far this has been proven only for fluctuations of 1/A given in the form of the gamma distribution, 

'(i)-i£y (£)"■»(-*)■ 

but this conjecture seems to be valid also in general J2Z] • 



3 Inelasticity obtained from analysis of the pp collider and pp 
fixed target data 

The available information in this case consists of: 
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(i) The mean multiplicity of charged secondaries, (n c h), produced in nonsingle diffractive reactions at 
given energy y/s, which can be parametrized as (n c h) = A + B In s + C In 2 s or as (n c h) = D + Es 1 
(see |28|V In what follows we shall assume, for simplicity, that they are pions with mass p = 0.14 
GeV. Out of it, we shall construct and use the total mean number of produced particles assuming 
it to be N = |(n c fc(s)). 

(ii) The observation that the phase space filled by the produced secondaries is essentially one-dimensional 
with limited (and only slowly growing with energy) transverse momenta |29| : (pt) = 0.3 + 

0.044 In ED]. 

(Hi) The rapidity distributions of charged secondaries, dN(s)/dy = f d 2 px^, provided either by col- 
lider experiments [141 1151 fTHl II Tj or by the fixed target experiment They are available only in 
a limited range of the rapidity space, depending on the details of the experimental setup |29| . 

Following the results obtained in J2j we expect (and therefore assume in what follows) that only a part 
W — Kyfs of the total energy y/s is used to produce secondaries in the central region of the investigated 
reaction. The inelasticity K = K(s) will therefore be the main quantity we shall investigate. 



We start with the information theory approach in its extensive version. In this case the relevant 
probability distribution defining information entropy is given by 



p(y) = p(y; N, W = Kyfs) = i • ^ 



(11) 



whereas constraint Q is just the energy conservation (here pr = \J P 2 + {pt) 2 and (E) is the mean 
energy per produced particle) |33| : 



dy [p T coshy] p(y) = = ( E) = K ^ ■ 



The limits of the relevant longitudinal phase space are 



Y m = In • 



w 

2p^ 



W' 2 J 



W = W — (N — 2)pq 



(12) 



(13) 



We would like to stress that throughout this paper the central region of the reaction, i.e., the region 
populated by produced particles distributed according to p(y) (or p q (y) later on), is always defined by 
cq. (|13|l . Therefore y £ (— Y m ,Y m ) in the CMS and we do not choose arbitrary cuts in rapidity space 
(as, for example, in |13jV Following now the steps mentioned in Section 2, i.e., minimazing the respective 
information entropy J3J) with the constraint given by (|12f) . we arrive at 



N dy 

with (3 — j3(W, N, pr) given by solving eq. fl2^l and 



1 dN(y) 1 
P\V) = — —J— = 7 exp(-^ T coshy) 



Z = Z(W, N, p t ) = J dyexp {-(3p T coshy) 



(14) 



(15) 



given by the normalization condition, J dyp(y) = 1 . Notice that in such an approach the " inverse tem- 
perature" (3 and the normalization Z depend on our input information, i.e., on N and pt, and on the 
assumed inelasticity K (via W = K ■ y/s), which is our free parameter. They arc therefore maximally 
correlated which means that the shape of distribution p(y) (given by (5) and its height (given by Z) are 
not independent of each other. This is in sharp contrast to the approaches presented before in [19] where 
both (3 (called "partition temperature") and the normalization (our Z) were treated as two independent 
parameters. Because of the symmetry of the colliding system momentum conservation does not impose 
any additional constraint here and we are left with f3 being the only Lagrange multiplier to be calculated 
from ea. 1)12(1 for each energy W and multiplicity N. 



4 



Notice that i|14|) . although formally resembling formulae obtained in thermal models [H21> has a much 
wider range of applicability as it is not connected with any assumption of thermal equilibrium. Actually 
it can be written in the scaling-like form: 



p(y) 



■ exp 



[It cosh y 

(E) 



(16) 



where (E) — W/N is the mean energy per produced particle and /? = /?• W/N. Plotting (3 as a function 
of (E) one observes that for the minimal number of produced secondaries (N — > 2) /3 — > — oo, whereas 
for the maximal number (N — ► N max = W/ht) [3 — > +oo. There is also an intermediate region in 
which /3 remains fairly constant leading to an approximate "plateau" in p{y). In this region "partition 
temperature" 1/(3 and inelasticity K are related in a very simple way, namely 







N 



K ■ 



(17) 



Actually, (3 > only for N > No = 2\n(N max ) and "plateau" occurs only for N ~ AT . It clearly 
shows that T = 1//3, called sometimes "partition temperature" ^Hjj is a measure of energy available per 
produced particle (which therefore depends also on inelasticity). 




Figure 1: Rapidity spectra obtained by UA5 
and Tevatron ^S] fitted by formula l|14|) . For com- 
pleteness results for pp data at 20 GeV [33 are also 
shown. The corresponding values of inelasticity K 
are listed in Tabled 



Figure 2: Rapidity spectra obtained by UA5 [TQ 
and Tevatron ^B] fitted by formula i|19|) ■ For com- 
pleteness results for pp data at 20 GeV [HJ are also 
shown. The corresponding values of q— inelasticity 
n q and nonextensivity parameter q are listed in 
Table [2 



In the case of the nonextensive version of information entropy the energy conservation constraint is 
given by 

£j dy [fiT coshy] ]p 9 (y)]o = ^ = *Lj^. (18) 
and maximization of the corresponding Tsallis entropy Q results in 

Pqiy) = JJ-^ = Y exp i (-^AT coshy) (19) 

(where exp g (...) is defined by eq. JHJ and Z q is given, as in (|15fl . by the normalization condition, 
J dyp q (y) = 1). The characteristic feature of p q (y), as shown in Fig. |SJ is that it enhances (depletes) 
the tails of the distribution for q > 1 (q < 1), respectively (or, in other words, it enhances or depletes 
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the more or, respectively, less probable events). Notice that in this case, differently than in Ijl4(l . one has 
to be sure that 1 — (1 — q)(5 q [iT coshy > 0, which imposes an additional condition on the allowed phase 
space for q < 1 and (3 q > 0. In this case the strict correlation between the shape of rapidity distribution 
and its height is relaxed because both depend also on the new parameter q. This fact will be important 
later on. The scaling- like formula (|16[) and the approximate relation (|17|l . this time between [3 q and n q , 
are still valid, albeit this time only approximately, i.e., for small values of \q — 1|. On the other hand 
P q > for N = (2h\N max ) q , i.e., for larger (smaller) multiplicities, depending whether q > 1 (q < 1). 



P238 data 
UA7 data 

K = 0.60, q = 1 .00 
K = 0.50, q = 1 .00 
k = 0.26, q = 1.27 
k = 0.2, q = 1.25 
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Figure 3: Comparison of data for rapidity distri- 
butions at a/s = 630 GeV obtained by P238 [TBJ 
(open circles) and UA7 1.7, collaborations (full cir- 
cles) with predictions of our model for it extensive 
(q = 1) and nonextensive versions. Notice that 
one can fit at the same time either P238 or UA7 
data but not both together. 



Figure 4: Energy dependences of inelasticities ob- 
tained in extensive (K, see l|12}l ) and nonextensive 
(K q , see ((2U) approaches (cf. also Table GJ. No- 
tice that results for P238 data do not follow the 
overall trend (see text). 



The results for fits using extensive formula l|14|l are displayed in Fig. Q]and those using its nonextensive 
version given by eq. I|19|) are shown in Fig. |2 As one can see they are almost identical, differences are 
showing up only for 20 GeV where data develop a tail which is most sensitive to the parameter q. The 
result for the joint distributions of UA7 data from ^Jj and P238 lfr data at 630 GeV are shown separately 
in Fig. Both sets of data are clearly incompatible in the sense that UA7 data (which are for 7r°'s 
and have been taken here in the same way as in |34| ) do not continue the trend shown by P238 data 
(which are, as UA5 ones, for charged particles). Instead they seem to continue the trend of the charged 
UA5 data, which is also clearly seen from the values of the obtained parameters displayed in Table ^ 
Because of this fact we have fitted them separately. The results for the respective inelasticity parameter 
K and its g-equivalent n q for different energies are shown in Tabled The estimated errors (the same for 
both approaches) range from AK = 0.02 for 20 GeV (where the fitted range in rapidity is biggest) to 
AK — 0.05 for 1800 GeV (where the lack of measured tails prevents a better fit). These errors should 
be kept in mind when looking at Fig. which summarizes the obtained inelasticities of different types. 
These inelasticities can be compared with inelasticities K\ y \<_^ defined by the formula 

N fYm=* 

K \y\<i = ~r dy p(y)[^ T cosh y} 7 (20) 

V s JY m =-i 

i.e., in the same way as in |19| . namely, by integrating over the same part of the phase space limited by 
|V m | = 4. Notice that values of K\ y \<± are systematically smaller than the corresponding values of K 
and are in a visible way decreasing with energy. The reason for such behaviour is that K\ y \ <4: counts the 
fraction of energy in a fixed domain in the phase space given by the condition that \y\ < 4 whereas our 
K gives the energy used for production of particles in the whole kinematically allowed region \y\ < Y m 
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as defined by eq. Table [I] contains also the corresponding values of the "partition temperature' 

T = 4 obtained from the eq. I|12|) . 
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0.45 


0.60 


36.29 


0.22 


1.27 


0.26 


35.51 
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0.22 
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31.0 


0.45 


0.50 


28.90 


0.21 


1.25 


0.20 


21.22 


0.40 


0.21 


900 


34.6 


0.48 


0.41 


29.47 


0.17 


1.29 


0.18 


30.79 


0.41 


0.17 


1800 


46.4 


0.50 


0.50 


55.69 


0.13 


1.33 


0.19 


62.57 


0.51 


0.13 



Table 1: The results for extensive (eqs. (11411 s ) and nonextensive (eq. (|19|| N 1 approaches applied to pp 
|31| and pp data |14l 1151 ITHl I17j on rapidity distributions. The first three columns summarize our input 
information: energy y/s, total charge multiplicities (n c h) (estimated as in [25]) and (pr) (estimated as 
in [Hj). Presented are extensive inelasticity K and parameters n q and q out of which the nonextensive 
inelasticity K q is calculated by means of (|21|) ). For completeness the corresponding "partition temper- 
atures" T — 1/(3 and r q — l//3 q are also listed as well as the corresponding inelasticities obtained for 
the limited portion of the phase space: ^|j,i<4 as given by eq. (EOJ and K q lvl - 4) as defined by (EH). For 
■s/s = 630 GeV we display separately results from fitting (a) P238 data ^| and (b) UA7 data |17j . 

In what concerns the nonextensive approach one must realize that parameter K q occurring in l|18|) 
is not the inelasticity in the same sense as K from the extensive approach, cf. eq. I|12|) . The reason 
is simple (and shown in best way in where Hagedorn statistical model of multiparticle production 
|3*6*] has been extended to (^-statistics) . Namely, because q summarizes all kinds of correlations and/or 
fluctuations present in the system (and make it nonextensive) the energies per particle present on the 
rhs of eqs. (|12fl and H18|) also contain in the nonextensive case contributions from these correlations or 
fluctuations, i.e., a kind of effective interaction characterized by \q — 1| The enhanced (depleted) 

for q > 1 (q < 1) tails of p q (y) observed in Fig. [5] say that more (less) particles are sent there towards 
the end of the phase space, respectively. This fact must then be compensated by the appropriate choice 
of the value of " g-inelasticity" parameter K q in (|18|l . which fixes the energy W q in this case. It should be 
also added at this point that when one is allowing the whole energy ^Js to be used for the production 
of secondaries and uses nonextensive version of information theory to find ^ then one finds q < 1 |37j. 
not q > 1 as here. This is, however, expected because of the already mentioned fact that the q < 1 
case enhances frequent events whereas q > 1 the rare ones. When considering the whole phase space it 
comprises all produced particles, which are located predominantly only in part of it. Therefore we have 
to enhance those frequent events by using q < 1. This choice, as was mentioned above, limits also the 
allowed phase space. On the other hand, when K is accounted for (as in the present case), the allowed 
phase space is essentially correctly described by K ■ -y/s and one has only to enhance the rare events when 
particles (because of fluctuations) "leak out" of it, what results in q > 1. 
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Figure 5: The examples of p q (y) as given by 
eq. l|19|) for q = 0.7 and q — 1.3 compared with p{y) 
for q = 1 as given by ea. ()14j) for one-dimensional 
hadronization of mass M =100 GeV into N = 20 
secondaries. 



Figure 6: The values of the nonextensivity pa- 
rameter q obtained in fits shown in Fig. |2] (and 
listed in Table Q} compared with the values of the 
parameter k of Negative Binomial distribution fit 
to the corresponding multiplicity distributions as 
given in |25|. see ea. lj23|) . 



In order to get a nonextensive version of inelasticity, i.e., K q , let us first observe that, according to 
<|12|) , inelasticity K can be expressed by the mean energy per particle, (E) , therefore, in the nonextensive 
case one can write accordingly 

N N f Ym 

K q = -7={E) q = —= / dy p q (y)[fi T coshy], (21) 

3-2q 

(with Y m provided by eq. (|13[l ). The approximate relation of K q with the parameters n q and q in 
nonextensive formula (|19f) arises when one estimates (E) q for \Y m \ — > oo and uses the nonextensive 
version of relation (|17|l . As it can be seen from Table nonextensive inelasticity K q defined this way 
agrees reasonably well with the extensive inelasticity K. Notice that K q v ^~ 4 \ defined as 

K (\y\<*) = dyp q {y)[n T coshy], (22) 

V s JY m = -4 

is essentially identical with K\ y \ <4: discussed above (and for the same reason). These results indicate 
that the true equivalent of inelasticity K in nonextensive approach is K qi at least in the sense used in 
cosmic ray research, namely that it defines the part of the initial energy taken away by leading particles: 
K elastic = 1 — K should be replaced by: K e i ast i C = 1 — K q . They are displayed explicitly in Fig. 0] Notice 
that P238 data a t y/s = 630 GeV clearly do not follow the trend presented by the UA5, UA7 and 
Tevatron data ^] ^j] . When neglecting this point the overall tendency is that inelasticity is essentially 
constant with energy and equal to K ~ 0.5, which agrees with first estimates made in j^Hl 01 1201 an d 
with first experimental estimates based on the analysis of the leading particle effect provided in |39| . 

A comment on the possible physical meaning of the parameter q obtained from our fits and listed 
in Table ^ is m order here. As we said before, in general, the nonextensivity parameter q summarizes 
the action of several factors, each of which leads to a deviation from the simple form of the extensive 
Boltzmann-Gibbs statistics, or Shannon entropy defined by eq.Q. Among them are the possible intrinsic 
fluctuations existing in the hadronizing system (23 ■ Notice that in the case considered here we have not 
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Figure 7: Fits to =p for different rapidity bins for i/s = 200 and 900 GeV (left and right panels, 
respectively) by means of ea. (fT9fl . 



accounted explicitly for the fact that each event has its own multiplicity N but we have used only its mean 
value, (N), as given by experiment where (N) = NP(N) with P(N) being the multiplicity distribution. 
Actually, we have used only its charged part, (rich), assuming that N — |(n c ^), i.e., neglecting in addition 
also possible fluctuations between the number of charged and neutral secondaries. Experimentally it is 
known that P(n c h) is adequately described by the so called Negative Binomial distribution (NBD) 29 , 
which depends on two parameters: the mean multiplicity {n c h) and the parameter k (k > 1) affecting its 
width, 

1 = cr 2 (n ch ) 1_ 

k (n ch ) 2 (n ch )' 

For k — > 1 NB approaches a geometrical distribution whereas for fc _1 — > it approaches a Poissonian 
distribution. In general it is found |29| that 

- = -0.104 +0.058 -In Vs. (24) 
k 

Following the ideas expressed in j2H] we would like to draw attention to the fact that the value of may 
be also understood as the measure of fluctuations of the mean multiplicity (cf. also |4()|). When there 
are only statistical fluctuations in the hadronizing system one should expect the Poissonian form of the 
corresponding multiplicity distributions. The existence of intrinsic (dynamical) fluctuations would mean 
that one allows the mean multiplicity h to fluctuate. In the case when such fluctuations are given by a 
gamma distribution with normalized variance D(n) then, as a result, one obtains the Negative Binomial 
multiplicity distribution with 



That is because in this case (see also |41p: 

h - 

r(fc) _ r(i + n)r(fc) (7 + i) fc +™ 



P „) = / dn— — ■ ^— — = — \ ■ (26) 
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where 7 = -Mr. Assuming now that these fluctuations contribute to nonextensivity defined by the 
parameter q, i.e., that D{n) = q — 1 [231 one should expect that 02] 

q = 1 + l- (27) 

As can be seen in Fig. this is precisely the case. Namely, fluctuations existing in experimental data 
for the rapidity distributions, dN/dy ^JESESL an d disclosed by fits using the nonextensive form <jl§|) 
follow (except for the lowest energy point at 20 GeV) the pattern of fluctuations seen in data for mul- 
tiplicity distributions P(N) and summarized by the parameter k of NBD (221- This means then that 
these data contain no more information than used here, namely existence of limited pr, inelasticity K 
and fluctuations as given by q > 1 or fc _1 > 0. 
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Table 2: Results for the same parameter q, K q , r q and K q as shown and defined in Tabled but now 
obtained for data taken in the restricted multiplicity bins AN ^B] and shown in Fig. Here: iVi = 76 
for Vs = 200 and 60 for 900 GeV. 

It is interesting to notice that whereas data on rapidity distributions ^] El H61 ITT] could be fitted 
both by the extensive (|14|) and nonextensive l|19|l distributions, similar data for rapidity distributions 
measured in restricted intervals of the multiplicity, AN |18| . can be fitted only by means of the nonex- 
tensive p q {y) as given by ea. l|19() 4cfJ. The extensive approach with maximally correlated shapes and 
heights of distributions, as discussed above, is clearly too restrictive. Only relaxing this correlation by 
introducing parameter q (i.e., by using nonextensive version of the information theory approach) allows 
for adequate fits to be performed, see Fig. [7|and Tabled Notice that partial inelasticities of all kinds are 
clearly correlated with the multiplicity bins, the higher the multiplicity the bigger is the corresponding 
inelasticity. The same kind of correlations are observed at 200 GeV between multiplicity and q, which 
increases with multiplicity. However, at 900 GeV q remains essentially constant. Following discussion 
in the previous paragraph one expects that it means an increase of the corresponding fluctuations. The 
question, however, remains, in which variable? We argue that the fluctuating variable in this case is 
inelasticity itself. The point is that particles filling a given interval of multiplicity AN can be produced 
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Figure 8: Inelasticity distributions x(K — K q ) (normalized to unity) obtained from ^2] data for s/s = 200 
GeV and y/s = 900 GeV (see text for details). K q is estimated from the q and n q parameters in Tablc|5]by 
using ea. (|2"T)l . To guide the eyes obtained results were fitted by simple gaussian and lorentzian formulas 
( X (K) ~ exp [-{K - (K))/(2a 2 )] and X (K) ~ a/[A(K - (K)) 2 + a 2 ], respectively, with a obtained for 
K e (0, 1)). Errors are estimated from the widths of the bins AiV and the shapes and errors of P(N) 
used in obtaining x(-ftT) 



in events with different, i.e., fluctuating values of K . As before, this fact would then lead to the apparent 
nonextensivity visualized by q > 1 and measuring also the strength of such fluctuations represented by 
the variance, 

a 2 (K) = (K) 2 (q-1). (28) 

However, in this case we do not have any independent estimation of &(K), therefore we could not exclude 
the action of some other, so far not yet disclosed, factors and propose the equivalent of eq. I|27|) for this 
case. On the other hand (|28(l could be used for estimation of the uncertainty in K once its mean value 
and the nonextensivity parameter are known. For example, taking from Table^the corresponding values 
of (K(s)) and q(s) one can estimate that (if, (200)) = 0.50 ± 0.23 and (if„(900)) = 0.41 ± 0.22. 

The results for partial inelasticities obtained from fits shown in Fig. allow us to calculate the cor- 
responding inelasticities K q (by using J3TJ), cf. Table These in turn, with the help of experimentally 
measured multiplicity distributions P(n c h) (taken in this case from |44| ) allow us to obtain, for the first 
time, the (normalized) inelasticity distribution x(K — K q ) presented in Fig. |H1 45 . This is one of 
the most important results which could be obtained only by using information theory approach in its 
nonextensive version. The gaussian and lorentzian fits shown here resemble very much the form of x{K) 
obtained in the so called Interacting Gluon Model of high energy processes developed and studied in 

(HUH urn 



4 SUMMARY AND CONCLUSIONS 

Using methods of information theory, both in its extensive and nonextensive versions, we have anal- 
ysed pp collider data ^] 1151 1161 1171 IT8"| and pp fixed target data [S] on multiparticle production. 
Our investigation was aimed at the phenomenological, maximally model independent description, which 
would eventually result in estimations of inelasticities for these reactions, their energy dependence and, 
whenether possible, also in the inelasticity distributions x(-^0- The information theory approach used 
by us leads to a comfortable situation where the only fitted parameter is either the inelasticity K (when 
using extensive approach) or parameters q and n q (in its nonextensive counterpart) out of which one can 
reconstruct the inelasticity K q or Kq~^\ It turned out, that data for rapidity distributions obtained for 
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the mean multiplicities can be fitted using both approaches. In this case the nonextensivity parameter 
q obtained from fitting the rapidity distributions dN/dy is practically identical to the parameter q, ac- 
cording to |23| , responsible for dynamical fluctuations existing in hadronizing systems and showing up in 
the characteristic Negative Binomial form of the measured multiplicity distributions P(N). This means 
therefore that in collider data for pp collisions and fixed target pp data analysed in this way, there is no 
additional information to that used here. On the other hand rapidity distributions for fixed multiplicity 
intervals AN can be described only by nonextensive approach and we argue that in this case q reflects 
fluctuations in the inelasticity itself. These data were therefore used to estimate, for the first time at this 
energies (200 and 900 GeV), the inelasticity distributions x{K), cf. Fig. |S] 

It is particularly interesting and worth of stressing here, that, formulas obtained by means of infor- 
mation theory are apparently identical with the corresponding equations of statistical models used to 
describe multiparticle production processes |32| . The point is, however, that - as was already stressed in 
appropriate cases before - the " partition temperature" T and the normalization Z are in our case not free 
parameters anymore. The only freedom (in our case it was the choice of inelasticity K) is in providing 
the corresponding constraint equations, which should summarize our knowledge about the reaction under 
consideration. Once they are fixed, the other quantities (in particular T) follow. This constrains seriously 
such approach and therefore in cases where it fails one can either add new constraints or include some 
interactions by changing the very definition of how to measure the available information. The Tsallis 
entropy used here [221 1231 1271 1261 13"5) is but only one example of what is possible, other definitions of 
information are also possible albeit not yet used in such circumstances [221 BE] ■ 
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